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Abstract 



O . Coherent Compton scattering on light nuclei in the A resonance region is 

studied in the impulse approximation and is shown to be a sensitive probe of 
the in-medium properties of the A resonance. The elementary amplitude on 
a single nucleon is calculated from the unitary K-matrix approach developed 
previously. Modifications of the properties of the A resonance due to the 
nuclear medium are accounted for through the self-energy operator of the A, 
\f} \ calculated from the one-pion loop. The dominant medium effects such as the 

Pauli blocking, mean-field modification of the nucleon and A masses, and 
particle-hole excitations in the pion propagator are consistently included in 
nuclear matter. 
2001 PACS numbers: 25.20.Dc, 13.60.Fz 
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2 1 I. INTRODUCTION 

P. ' 

Coherent Compton scattering on nuclei in the region of the A resonance is of considerable 
interest. The reaction allows one to study the propagation and decay of the A in the nuclear 
medium. In particular the shift of the pole position and a change of the width of the A 
peak, reflect sensitively in the cross section and polarization observables. 

For a comprehensive review on Compton scattering we refer to the recent reference [0. 
Here we only mention several models relevant for the present study. The so-called " schematic 
and extended" models 0j3| made use of the optical theorem and dispersion relations to relate 
the Compton-scattering amplitude with the total photoabsorption cross section. Important 
contributions to the amplitude, such as nuclear kinetic and meson seagull terms, and A 
excitation, were added explicitly ||. The data for Compton scattering were described quite 
well at forward scattering angles. 

Compton scattering was extensively studied in the A-hole model |3r|| which was origi- 
nally developed by Koch, Moniz and Ohtsuka in MM. The A excitation is included through 
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the Ml transition, the background is represented by the Kroll-Ruderman part of the virtual- 
pion photoproduction amplitude, while in Refs. IHH the proton Thomson term (with El 
multipolarity) is also added. 

The starting point of the present model is the full Compton-scattering amplitude T N 
on a free nucleon that includes all dominant mechanisms at energies up to, at least, the A 
resonance. This amplitude has been obtained in the unitary K-matrix approach, similar to 
the calculations presented in refs. [|5]|| for pion-nucleon scattering, pion photoproduction 
and Compton scattering on the proton. Parameters of the model are fixed from a detailed 
fit to observables for the above elementary reactions. The Thomson limit dictated by the 
low-energy theorem on the nucleon and the Kroll-Ruderman term of the photoproduction 
amplitude are automatically included. 

Medium effects are taken into account via the self-energy of the A, while the self-energy is 
calculated from the pion-loop mechanism embedded in nuclear matter, as described in detail 
in sect. II. In particular, Pauli blocking effects, mean- field modifications of the nucleon and 
A masses, and particle-hole contributions to the pion propagator are consistently calculated. 
Thus we avoid arbitrariness associated with introducing a parameterized self-energy. 

The amplitude on a nucleus is calculated in the impulse approximation. A formfactor is 
introduced to account for finite size effects. The specific density dependence of the medium 
effects is accounted for through a local density approximation. Fermi-motion is treated 
in the so-called factorization approximation, in which the elementary Compton-scattering 
amplitude is calculated on a nucleon with an effective momentum [p!0| -[T3|. The latter is 
taken in such a way that the energy-momentum conservation for the 7 N scattering holds, 
thus ensuring gauge invariance at the one-body level. 

Differential cross sections and photon asymmetries are calculated for light nuclei ( 4 He 
and 12 C) in the energy region 150 - 250 MeV. 

The structure of the paper is as follows. In sect. II we describe the calculation of 
medium effects in the A propagator. Sect. Ill presents the basic formalism for calculat- 
ing the Compton-scattering amplitude on nuclei in the impulse approximation. Results for 
the cross section and photon asymmetry are presented in sect. IV. Conclusions and prospects 
are further outlined. In Appendix A formulas for the pion self-energy are collected. Details 
of the kinematics involved in the reaction on a nucleus are contained in Appendix B. 

II. MEDIUM MODIFICATIONS OF THE A RESONANCE 

The properties of the A in the nuclear medium are calculated in a relativistic framework 
for symmetrical nuclear matter, along the lines of Refs. |fH|-|T(][] . The medium modifications, 
which are expressed through the dressing of the A propagator, are investigated using different 
levels of approximation. The imaginary part of the A self-energy (or the A decay width) 
is calculated in different models for the nuclear medium. Dispersion relations are used to 
determine the real part (mass modification) of the A self-energy. 

In a first subsection the Rarita-Schwinger formalism is introduced to describe the A 
resonance. The next three subsections deal with the A decay width. In a first step the nuclear 
medium is presented as a non-interacting Fermi-sea of nucleons leading to the calculation of 
the Pauli-corrected decay width and the spreading width. Subsequently mean-field effects 
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are included within the au— model of Walecka JT7j. Throughout, the conventions of Bjorken 
and Drell [18] are used, and all expressions are derived in the nuclear matter rest frame. 



A. Formalism for the A self-energy 

The free A propagator in the Rarita-Schwinger formalism has the following form [I9|: 

(1) 



G A (p A ] 



^ A - M A + it 



3 ' ' 3M A 3M A 



with the mass Ma = 1232 MeV and p A = (p A , p A ). The invariant mass of the A resonance 
will be denoted by W A = ip 2 A ) 1 ^ 2 - 

The propagator has to be dressed due to the interaction of the A with the nuclear 
medium. This is taken into account by introducing the self-energy in the equation for the 
inverse propagator 

(G^r = ((Giy'r - = - m a ) (v 3/2 r 

- 2(j* A - M A ) (Pn /2 r + ^M A [ {V^T + {V^iT (2) 

where the inverse free propagator is written in terms of the projection operators on the spin- 
3/2 and spin- 1/2 sectors [|(J (P 3//2 , V]j 2 , ij = 11, 22, 12, 21 ). In this paper we concentrate 
on the dominant spin-3/2 component of the A propagator, and therefore the self-energy is 
taken as 

E A " = T?l 2 (V ZI2 Y V , (3) 

where we discarded terms proportional to the spin-1/2 projection operators. The spin-3/2 
component of the self-energy has the structure 

S a /2 = ^a(pa) + j Da(pa)^a, (4) 
and the dressed propagator becomes 
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W = r— - 375 r /2 ) 

p A -M A - 



= [^a(1 - ReD A (p A )) - (M A + ReC A (p A )) 

- i(ImC A (p A ) + ImD^p^)}- 1 (V 3/2 r. (5) 

Note that in vacuum C A (p A ) and D A (p A ) are functions of the invariant W A only, while in 
the medium these functions acquire additional dependencies on the 3-momentum \p A \ and 
the nuclear density p. We keep p A as an argument in Eq. @ for the general situation. 



B. Width of the A resonance and Pauli blocking 

The imaginary part of the A self-energy is calculated from an effective interaction La- 
grangian. We focus on the pion-nucleon decay channel in a nuclear environment. The 7riVA 
interaction Lagrangian is taken in the form |20|-|22l 
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C.NA = — g 1 "^* {d v w) + H.c. , (6) 
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where \1/ (A M ) is the nucleoli (A) field operator, n the pion field, and % is the isospin | 
transition operator with the normalization condition TjT^ = 5y — , (i, j = 1, 2, 3) . In 
the above Lagrangian we dropped the off-shell coupling of the A described by the parameter 
z v |]20|-p2|, since it affects only the spin-1/2 components which are not considered in the 



present work. 

The width of the resonance is defined by the imaginary part of the self-energy as 

T A (p A ) = -2Im[C A ( PA ) + W A D A (p A )}. (7) 
This can more formally be written trace 



Ta(pa) = -21 m 



4 ir LA (? W 2W A 



(8) 



The contribution to the self-energy for a A decaying into a nucleon and a pion is 

'/ttWAX 2 ^ f d 4 k 



*(£a)^=(^) 2 E I ^%TlG N {k){-q v )T a DM, (9) 



a=l,2,3 ' 



where q = PA~k, and Gjy(k) and D n (q) are the nucleon and pion propagators, the structure 
of which depends on the model used for the nuclear medium. We will now discuss different 
approximations and their implications for the A decay width. 

In vacuum the nucleon and pion propagators in @ reduce to the free nucleon propagator 

G° (k) = 1 = Me ( A+ Ck) _ A-(-fc) \ 

N{ } fc-M N + ie E % [k -E % + ie k Q + E % -ie) 1 ] 

and the free pion propagator 

D°M = - ^— , (11) 



where E% = y\k\ 2 + and A ± are the projection operators 

Applying the Cutcosky rules for the discontinuity of the loop integral in Eq.flj^), and 
using 2o=i 2.3 TlT a = 1 one gets the conventional expression for the width in vacuum 



fn NA (Ek, + M N) , , 
127rm2^ A 



-ptiac _ JnNA V-^/Ctt 1 _^J_ 1^ (13) 

where 
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4w A 



[Wl - {M N + m^f}[Wl - (M N - m n ) 2 } (14) 



and E^ = y ^1 + ^n- Putting = Ma in Eq. ( p~3| ) and using the experimental value 
pwtc — 115 M e v one obtains the coupling constant f W NA = 2.15. 

The effect of Pauli blocking is included by replacing the free nucleon propagator by a 
medium- modified one. If one models the nuclear environment as a non-interacting Fermi-sea 
of nucleons the propagator becomes (see, e.g., 



G F N G W 



\k\-k F ) 6(k F -\k\) 




M N - ie 



8(k°) + 




W-*,) + ^WI,_^^, (15) 



fco — + ie k — E% - ie 
This is written in the nuclear-matter rest frame, an explicitly covariant form is discussed 



m 



14]. Implementing this propagator in Eq. (0) and defining the Fermi-energy Ep 



k F + Mfj and E± = (p° A E k7T ± \pA.\k n ) /W& we obtain the Pauli-corrected A decay width 



^a(Pa) = J *"%\ **' I * [E + - max(£_, E F )\ 6(E + - Ep) . (16) 

This result can be separated into three different energy regions: 

1. E_ > E F : the Pauli-corrected decay width reduces to the vacuum expression 
making the width (-independent; 



2. E,>E F > 



^ = + WAlk ' ~ W * Ef) • (17) 



The width becomes \pa\- and density-dependent. Taking the limit \p&\ — > reduces 
the width to its vacuum value if E^ > Ep, and to zero if E^ < Ep] 

3. Ep > E + : the width becomes zero (Pa(Pa) — 0). This means the Fermi-sphere 
engulfs the decay sphere completely, making the Pauli blocking complete. 

In Fig. (0)(a) the full dependence on p and \p&\ is shown for = 1232 MeV. The results 
at densities 1.2, 0.8, 0.4 and 0.05 times normal nuclear matter density po (kp — 1-333/m -1 ) 
are plotted separately in Fig. (|l])(b). At high A-momenta and low densities the energy of 
the decay nucleon lies well above the Fermi energy, and no blocking occurs. At somewhat 
lower momenta of the A part of the momenta of the decay nucleon are Pauli blocked. With 
increasing density this blocking may become complete for the lowest A-momenta making 
the A unable to decay into a pion-nucleon pair (see dashed and dotted curves in Fig. (|l])(b)). 
These phase-space considerations result in a strong energy-dependence of the A decay width. 
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C. Spreading width 



In the nuclear medium the pion will strongly interact with the surrounding baryons 
creating nucleon-hole and A-hole excitations. This can be taken into account by dressing 
the pion propagator with the proper pion self-energy 



DM 



1 



q 2 - ml- U*(q) + it 



(18) 



where 11^ (g) = U. p h(q) +H-Ah(q) is the polarization self-energy of the pion. In our calculations 
we use the pion-nucleon pseudo- vector coupling 



C 



nNN 



NN 



(19) 



with the 7riV./v"-coupling constant f n NN = 1-01 [13]. In this work we limit ourselves to 
forward and backward scattered particle-hole excitations, and omit anti-nucleon excitations 
in Eq. (|T5|), which play a role only for very large pion momenta. Also intermediate A- 
hole states are omitted from the pion self-energy. We expect that their contribution is 



more suppressed than the estimate in ref. |14| when the width of the A resonance is taken 
into account. In principle, a complete calculation of the A-hole states would require self- 
consistency between the pion and A self-energies, which falls outside the scope of the present 
paper. 

The lowest-order pion self-energy now reads 



n° 



ph 



(?) 



NN 



AE^ E^ 

^k+q^k 



d A k m 

Tr 

(2vr) 4 

0(\k + q\ 



— kp) 9{k F -\k\) 



+ 



(fc + Qo - + «e) (k -E % - 
9(k F -\k + q\) 6(\k\-k F ) 



it 



E t 



?e 



NN 



(k + q ~ - ie) (k 
\q\ 2 C (\q\)~C 2 (\q\) + 2M 2 N q 2 L (q) 



(20) 
(21) 



where k = (E%,k) and I = (E% + g,k + q). The analytical expressions for the functions L , 
Cq and C*2 are given in Appendix A. When summing the series of particle-hole bubbles 
the effects of short-range correlations are important. These short-range correlations are 
accounted for in the standard way by introducing the Landau-Migdal parameter g' NN = 0.6 



24 



n P h(g) = q 2 



njfc(?) 



q 2 + g' NN U° ph (q) ' 



(22) 



Using the pion self-energy IL^ in the pion propagator from (18|) we get the spreading width 
of the A resonance in the medium 
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^ - (m) fe)7^^ft£^(*-^^: 

(w 2 a 2 U n\A 1 ImU ph (q) 6(-q 2 ) 

[W A q - {pa ■ q) ) — f/}2 _ 2 _ „ n -™ - r . n -™ (23) 



A, 

X 



2^ [<f -mj- i?ell p/l (g)] 2 + [JmlV 



where q = p& — k. Since the pion self-energy only receives contributions from the particle- 
hole excitations, its imaginary part is non-zero only for space-like pion momenta. This puts 
restrictions on the integration boundaries in Eq. fl2~3"|). The two-dimensional integral was 
evaluated numerically. As an example we have plotted in the upper panel of Fig. (0) the 
imaginary part of the pion self-energy in function of the pion momentum |g|, as it appears in 
Eq. ( f2~3|) for the specific A kinematics : = 1232 MeV, \p&\ = 200 MeV and pion energy 
q° = (p% - E F )/2. 

The results for the spreading width are shown in Fig. @. The spreading width is roughly 
proportional to the density, which can be understood on the basis of the phase space available 
for the hole states. As can be seen from Fig. @(b) it is only weakly dependent on the 3- 
momentum | Pa |- Also the dependence on Wa turns out to be rather weak in the region 
of interest. The total width of the A in this non- interacting Fermi-sea of nucleons is given 
by the sum of this spreading width and the Pauli-corrected decay width from the previous 
section. 



D. Mean-field effects in the nucleon and A self-energy 

A refinement to the free Fermi-gas model can be made using the Walecka <juj— model |T7 



in the mean-field approximation. Here the a— and uo— meson couple to the nucleon resulting 
in the mean scalar and vector fields ($ s ) and (V^). 

In the nuclear matter rest frame the spatial part of (V 1 ) is averaged to zero, and the 
constant mean-field contribution to the nucleon self-energy becomes 



M*O = -&S + 0«^7o, (24) 

where and are the coupling constants of the scalar and vector field with masses m s 
and m v respectively, and ($ s ) = $, (V^) = 5^qV. This self-energy can be implemented in 



the modified Dirac equation [[17]] 



[]k-M N -E N (k)]V N (k) = 0. (25) 

Introducing the effective nucleon four-momentum k* = (k° — g^V,k) and mass = 
Mn — g^&, the nucleon spectrum is modified to k* 2 = M^ 2 , or ko = g^V±E*>, where 



Ei = ^\k\ 2 + (M N 



2 + (M* N ) 2 . 

In order to assess the sensitivity of the results to the mean-field parameters we have 



performed calculations taking 2 parameter sets from [23], henceforth called set I and II. 



Set I, called QHD-I in ||23|| , results from a pure mean- field approximation to the bind- 
ing energy. The dimensionless ratios of coupling constants and meson masses have values 
C 2 = (g^M N /m s ) 2 = 267.1, C 2 = (g*M N /m v ) 2 = 195.9. The nuclear matter equilibrium 
density is at k° F = 1.42/m -1 , with binding energy 15.75 MeV and an effective nucleon mass 
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M^/Miy = 0.56 at po- The full density dependence of the effective mass is defined explicitly 
by the self-consistency equation 



M* N = M, 



N 



N 



nit 



Ti 



k F E F -{M* N ) 2 \n 



k F + Ei 



M 



N 



(26) 



k% + (M* ) 



Set II, called the relativistic Hartree approximation in 



where E F = 

takes into account vacuum fluctuation corrections to the binding energy. The parameters 
are C 2 = 228.2, C 2 = 147.8. The equilibrium density is taken at k F = 1.30/m -1 , with 
a binding energy of 15.75 MeV leading to an effective nucleon mass M^/Mn = 0.73 at 
equilibrium density and the self-consistency equation reads : 



M* N = M, 



N 



N 



+ 



(9 



ml 

N\2 



Ti 



k F Ep 



M 



N 



ml rr 2 



M* N 6 \n 



M 



N 



M, 



N 



- M 2 N (M* N - M N ) 
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- -M N (M* N - Mjv) - -(.!/; - M 
2 6 



N 



(27) 



The full density dependence of the effective nucleon masses in both cases are shown 
in Fig. (|^). We see a strong reduction of the effective nucleon mass compared to the free 
nucleon mass with increasing density in both cases, with a slower decline when the influence 
of negative energy-states is considered. 



In the extended mean-field model of ref. [141 the A is assumed to move in the mean a and 



ijj fields. The mean-field contributions to the A self-energy can be treated in an analogous 
way as for the nucleon, i.e. they are absorbed in the effective A mass and 4-momentum 
p\. In this paper we employ the so-called universal couplings |TJ]: = gf- and g^ = g^ 
and as a result the A effective mass M£(p) may be expressed as M\{p) = Ma — (Mjv — M^). 

We can now investigate the influence of these mean-field modifications on the Pauli- 
corrected decay width and the spreading width of the A by replacing the Fermi-gas nucleon 
propagator (Eq. (|T5|) ) in the expression @ and in the calculation of the pion self-energy 
(see Eq.(|20D) with the mean-field propagator 



M^ + ie 



+ 



0(A*-|A;|) 



M 



N 



- + ie 



M 



N 



Et 



A* + (k) 



K\k\-k F ) 
h£-E$ + ie 

v k 



+ 



Q(k F -\k\) 



-k) 



kt-Ei 



ie 



kt + Ei- ie 

u k 



(28) 



where A* 1 * 1 follows from Eq. flT2] ) with the nucleon energy E^ and the nucleon mass M/v 
replaced by the effective variables Ei and Mjy-. The resulting expressions are formally 
identical to the free Fermi-gas expressions, if all kinematical quantities are replaced by their 
effective equivalents. In particular, we introduce the effective in-medium mass of the A, 
Wl = {p\ 2 ) l/2 . In what follows we will use these effective kinematical quantities in all 
expressions; it is understood that they reduce to the original kinematical quantities for the 
free Fermi-gas calculations. 
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In the middle and lower panel of Fig. (Q) we depict the imaginary part of the pion 
self-energy in function of the pion momentum |<f| for both mean- field calculations at the A 
kinematics : W A = M A , \pa\ = 200 MeV and pion energy q° = [(p A )° — E F \/2. 

The results in the era;— model for the Pauli-corrected decay width and the spread- 
ing width F A at the on-shell point W A = M A (p) for both parameter sets are depicted in 
Fig. (H) and Fig. (|6|). It is seen that the structure of the decay width hardly changes when 
effective masses are introduced; only the limiting value at large | p A | now becomes density- 
dependent. Because of the stronger reduction of the effective masses the Pauli-blocking is 
more pronounced using parameter set I. 

The mean-field effects result in an overall reduction of the spreading width as compared 
to the Fermi-gas calculation (see Fig. (j^)). For the relevant nuclear densities p/po < 1.2, this 
reduction is stronger at larger densities. It can be shown that for vanishing A momentum 
the spreading width is roughly proportional to the integrated pion propagator divided by the 
effective A mass. The reduction of the effective mass explains the global density dependence 
of the spreading width in this region. For larger densities the effect of the pion propagator 
makes the spreading width saturate and eventually decrease in the mean-field models. The 
mean- field model I yields a maximal spreading width at around the equilibrium density. In 
the mean-field model II, the spreading width saturates at much larger densities. 



E. Real part of the A self-energy. Renormalization 

In the previous sections we have obtained the imaginary part of the A self-energy due to 
the pion dynamics, i.e., the sum of the Pauli-corrected decay width and the spreading width. 
This imaginary part generates a contribution to the real part of the A self-energy, which can 
in general be obtained via a dispersion relation. Based on the general structure of the A self- 
energy in Eq. (H), one can find Re[C A {p A )] and Re[D A (p A )] through a dispersion relation 
in terms of Im[C A (p A )} an d Ifn[D A (p A )} calculated for arbitrary W A . The imaginary parts 
of C A (p* A ) an d D A (p* A ) are retrieved from the following relations 



Im C A {p\) 



Im D A {p* A ) 



AW A 



(29) 



where T A is defined as 



1 A 



-21m 




(p3/2) 



(±f A + wt) 



2W\ 



-2 Im[C A (p A )±W A _D A (p* A )}. 



(30) 



For instance, the imaginary parts for the one-pion loop in vacuum (W A reduces to W A ) 



arc 



Im C A (p A ) 



N 



2AnmlW A 



Im -D A (pa) 



2AnmlWl ' 



(31) 



The in- medium expressions for ImC A and ImD A (at non-zero nucleon density) are more 
complicated than Eqs. (|31~D , and depend on both W A and \p A \- 
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We make the assumption that an unsubtracted dispersion relation holds at fixed values 
of \pa\, namely 

Re{C A (W^, |p a |) + f A D A (Wf, |p a |)] = 

1 [to CJW k \ \ M + Kto D^Wl?, IP.I)] , 2 2 

where the formfactor /a(Wa 2 ) i s introduced for convergence. Eq.([3^) actually implies two 
separate relations for Ca and -Da- The threshold W*/, is M^r + in vacuum, while in 
medium it is a more complicated function of masses, Fermi momentum, and 3-momentum 

The propagator is renormalized in such a way that in vacuum it has a pole at the physical 
mass Ma = 1232 MeV with a residue equal to unity, as if the A were a stable particle. These 
conditions give rise to the renormalized functions (for arbitrary Wa) 

ReCl{pA) = ReC A {p A ) - SM A + M A (1 - Z' 1 ) , (33) 
ReDl{pA) =ReD A (p A )-(l- Z- 1 ). (34) 

The mass shift (the physical mass minus the bare mass) and wave function renormalization 
constant are given respectively by 

5M A = [ReC A (p A ) + M A ReD A (p A )]\ WA=MA , 

Z = {1- ReD A (p A )\ WA=MA - 2M A [ReC' A {p A ) + M A ReD' A (p A ) ] | Wa =m a } _1 , (35) 

with the notation ReC A (p A ) = dReC A (p A )/dW A and similarly for ReD' A (p A ). The required 
properties of the renormalized propagator in vacuum are ensured by the relations 

^/2\R\ = d , 3/2-.R, 

/a=Ma <9^A '/a— Wa 



(S 3 a /2 ) H = C|(p a ) + Cf(p a )^ a . 



Finally, the in-medium A propagator, which is used in the calculations described in 
sect. |TJ, reads 



GX(Pa) = [A (l - ^a(Pa)) - ( m a + i2eC2(?4) 

(/mC|(pl) + JmD2(pi)^)] _1 (V^ 2 r ■ (37) 



—2 



In the calculation of the dispersion integrals we use a similar formfactor as in [25 



f lw > 2 i - ( ^i-K) 2 + (Mi-M},r \ 2 

f - l(Ai-A W + (MJ-H-l)^ (38) 

where W A = (W A + M A — M A ) 2 and the normalization is chosen such that f A (M A ) = 1 
which is appropriate for the in-medium calculation. The cut-off parameter is taken the same 
as in 



25j: Aa = 1.506 GeV. 
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The mean-field description of the nuclear medium necessitated another modification of 
the formfactor. At large values of the nuclear density the decrease in the effective masses 
of N and A results in values close to zero for the threshold invariant mass W t h of the 
dispersion integral (|32|) , if \p A \ is large. Since the projection operator on the spin-3/2 sector 
of the A appearing in the calculations of r A , leading to Im[C a(pa)] an d I™[D&{p\)}, nas 
negative powers of the A invariant mass, this would cause unphysical large contributions 
to the dispersion integral coming from the region of the invariant mass close to zero. We 
eliminated these contributions by multiplying the formfactor in Eq. fl38| ) with g(W A ) = 
8(M N - W A )(1 - (W A - M N ) 2 /M* N 2 ). We checked that the multiplying factor g(W A ) hardly 
changes the A real self-energy for the values of \p A \ and p that enter the description of the 
Compton cross section in the next section; it is added simply to extend the p and \p A \ range 
of validity of the real part of the A self-energy. 



III. COMPTON SCATTERING 

The amplitude for the process of Compton scattering on a finite nucleus is calculated in 



the impulse approximation. We apply the so-called factorization approximation (see |26 
ch.ll, sect. 2) which was shown to work well in pion photoproduction |iOHl2| and pion 
scattering [0^7| on nuclei, in particular for the light nuclei in section IV, where the nuclear 
wave function is well described by a harmonic oscillator model. A large part of the effects 
of the Fermi-motion are accounted for by evaluating the amplitude on a nucleon moving 
with the effective momentum p (p f = p + q) in the initial (final) state, where q = k — k' 
is the momentum transfer (see Appendix B for precise definitions). The amplitude in this 
approximation is written as 

K A = A(T N (p))F p (q), (39) 

where F p (q) is the Fourier-transform of the density distribution (formfactor). In Eq. fl39|), 
the formfactor of the ls-lp-shell nuclei with Z = N = A/2 is constructed on the basis of 
the experimental charge densities in |28| (see table V therein), correcting for proton finite 
size effects and assuming equal proton and neutron densities. (TV) is the spin averaged 
single-nucleon amplitude defined as 

(T N (p)) = - <k\\'^fi,m s \T N \k,\^p,m s >, (40) 

m s =±l/2 

where p is the effective nucleon momentum in Eq. flBH| ) and m s is the projection of the 
nucleon spin on the OZ axis. For isospin-saturated systems (to which we restrict our present 
discussion) an isospin average is also performed. 

The cross section for unpolarized photons in the Lab frame can now be expressed as 

do_\ k ' 3E A 1 v \ K ,2 uu 

cM' lLab ~ 4Tr*kM A 2^, 1 M ' 1 ' 

where E' A = u + M A , uj = k — k' and A 7 , are the photon helicities. It is convenient to 
redefine the Compton T- matrix through the amplitude T^{p) 
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/ M N M N 1 lxi/2 
rjv(p1 = U jy (p)^(p>) 2fc2^) ' (42) 



The latter is normalized according to |29[ (App.A.3) and has simpler properties under 
Lorentz boosts. Calculating the amplitude in the jN CM frame (marked with superscript 
"c") we obtain 



A 7 ,A 7 



The photon asymmetry which can be measured with a linearly-polarized photon beam 
is defined as 

£ 7 = d ^~ d /W , (44) 
7 d a ± + d a\\ y ' 

where do\\ (da±) is the cross section for the photon polarization vector in the scattering 
plane (perpendicular to it). 

The cross section in the 7 A center-of-mass frame can be related to the cross section in 
the Lab frame as 



T7 -^L- (45) 



da _ / k \ 2 M\ da 
lcm= \~k>) 77 dfi' 

where sa = Ma{Ma + 2k) is the total invariant energy squared, and the center-off- mass 
photon momentum and the scattering angle are 

k cm = kM A /^/sA, (cos6» 7 ) cm = 1 - — (1 - cos6> 7 ) . (46) 

"'cm 

The single-nucleon amplitude is decomposed into one part which corresponds to the 
amplitude on the free nucleon, plus a term which accounts for the modification of the A 
resonance in the medium, i.e. 

T N = T } N ree + (Kfr - Kp) . (47) 

The first term is the T matrix for Compton scattering on the free nucleon; the term between 
brackets accounts for the nuclear-medium modification of the A resonance. To avoid double 
counting the vacuum contribution is subtracted. 

The T-matrix for Compton scattering off a free proton, T^ ree , is calculated in a K- 
matrix model very similar to that of ref. ||. This covariant coupled-channels calculation 
of pion scattering, pion photoproduction and Compton scattering on the nucleon, satisfies 
unitarity constraints below the two-pion production threshold and is gauge invariant. In 
the calculation of T^ ree the A is treated as a genuine spin-3/2 resonance |30| in order to 



be compatible with the present treatment of the in-medium A resonance. The change in 
the structure of the 7./VA and 7riVA vertices (the disappearance of the spin-1/2 off-shell 
couplings) necessitated modification of parameters of the p and a exchanges in the t - 
channel. A comparable fit to the data as in ref. || could be obtained. In Fig. (|7]) the 
results for Compton scattering are compared to data. At the pion-production threshold, 
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E 7 fa 150MeV, the calculation overestimates the data which might be related to ignoring in 
the K-matrix calculation the real pion-loop contributions which are responsible for the cusp 
structure in the f]f E Compton multipole p2| . |33| . 

The dressed A contribution K^ d is based on a calculation in which only the s-type tree- 
level contribution is taken into account using the medium-modified A propagator as defined 
in sect. II (see Eq. (|37|)). Note that in the impulse approximation the photon is absorbed 
on a free nucleon and thus one has to work with the free nucleon mass instead of the 
medium- modified one M^. The A self-energy parameters will depend on the difference 
between the A invariant mass and the nucleon mass (5 = — M^) and are therefore 
evaluated at = Mjy + 5. The subtracted A contribution K N } in vacuum is obtained 
from a similar calculation using the free propagator instead. 

In the limit of low photon energies the cross section for Compton scattering is given 
by the Thomson limit, where the matrix element is proportional to Z 2 / A. In the present 
calculation only the contribution to Compton scattering proportional to Z, the total number 
of protons, is taken into account, thus the contribution proportional to NZ is omitted. As 
such the Thomson limit is violated since Z 2 j A = Z — NZ/ A. The neglected contribution 
is thought to arise from intermediate excitations to collective giant dipole resonance states 



and be related to the finite extent of the nuclear system p4| , p5| . For this reason one expects 



this contribution to be vanishingly small at forward angles. At backward angles, where 
the one-proton cross section is suppressed by the formfactor due to the large momentum 
transfer, the two-body mechanism may give a significant contribution. In our approach 
the equivalent contribution would arise from two-body contributions to the electromagnetic 
current arising from the nucleon-nucleon interaction. As argued, such a contribution is of 
marginal importance at forward angles but large at backward angles. In a future work this 
will be included explicitly, presently we have ignored these two-body currents. 



IV. RESULTS FOR COHERENT COMPTON SCATTERING 

Cross sections have been calculated for 4 He and 12 C at several densities to investigate 
medium effects. To compare with data an average over density (pa), based on the Local 
Density Approximation (LDA), has been performed. The density profile (pa) was taken 
consistently with the formfactor in Eq. (|39|) . 

In Fig. (||) we have plotted, for various nuclear densities, the cross section and photon 
asymmetry for Compton scattering on 4 He in mean-field model I, both at fixed lab angle 
6^ = 37° and lab energy E 1 = 206 MeV. The results show a strong density dependence. In 
order to obtain more insight we have plotted in the upper panels of Fig. ([|) the values of the 
3- momentum | and (kinematical) invariant mass of the A as enter in the calculations 
presented in Fig. (|8]). In the lower panels we show the real and imaginary part of the A 
self-energy. We concentrate on the dominant imaginary part, as this seems sufficient to 
explain the global density dependence of the cross section. At fixed E 1 = 206 MeV by far 
the largest contribution to the imaginary part is due to the spreading width. The decay 
width in this energy regime vanishes for the larger densities and is very small for the lowest 
densities. The width therefore almost vanishes at zero density. At fixed 9 1 = 37° and with 
increasing energy the Pauli-corrected decay width becomes more important, showing up in 
the A width at small density and in the global increase with energy starting at 300 MeV. 
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Much of the density dependence of the cross sections in Fig. (H) can be understood from 
the density dependence of the imaginary part of the A self-energy. At a photon energy of 206 
MeV one is relatively far from the peak of the A resonance. An increase in the width of the 
resonance therefore results in an increase of the cross section at this energy. The opposite 
happens when one approaches the peak of the resonance, where the cross section decreases 
with density. The data show clear evidence that this is indeed the correct mechanism, at 
206 MeV the vacuum calculation falls below the data while the LDA result shows a good 
correspondence with the data at forward angles. Near the resonance the vacuum calculation 
overestimates the data by a factor 2 while the LDA result gives a much better prediction 
or even lies below. The sharp fall-off of the cross section with angle is mostly due to the 
formfactor which falls off strongly with increasing momentum transfer. 

At backward angles the cross section is not reproduced, which is probably due to the 
double-scattering contribution which is missing from the present calculations. The photon 
asymmetry at 206 MeV shows only a minor density dependence as compared to the error 
bars on the data. 

In Fig. ( |T0"D we compare the 4 He cross section and asymmetry with LDA calculations for 
the Fermi-gas and mean-field calculations I and II. The Fermi-gas calculation undershoots 
the data at small angles for E 1 = 206 MeV and at large energies for # 7 = 37°, and deviates 
from the asymmetry data points. The mean field calculations tend to improve this. 



The cross sections for C is shown in Fig. (11). Because of the larger radius of C the 



cross section falls off faster with angle than that for 4 He. The drop in the cross section at 
energies beyond 250 MeV is partly due to an increased width of the A resonance and partly 
due to the formfactor cutting the cross section at larger momentum transfers. This effect is 
also seen in the data. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have presented a calculation of the cross section for Compton scattering 
on 4 He and 12 C in a K-matrix model where the amplitudes are calculated in the impulse 
approximation. Fermi-motion is incorporated using the factorization approximation scheme. 
The medium effects are included by replacing the free A propagator by a medium-modified 
A propagator in the s-type resonant tree-level diagram. The medium-properties of the A 
resonance are investigated in a relativistic framework for symmetrical nuclear matter. 

This involves the calculation of the A self-energy which was performed in a Fermi-gas 
model and two mean-field models. The imaginary part includes the Pauli-corrected decay 
width and the spreading width incorporating ph-excitations; the real part was calculated 
using dispersion integrals. In both the Fermi-gas and mean-field calculations the width 
is increased as compared to the free A width. This increase tends to be stronger for the 
Fermi-gas model than for the mean-field models. 

The differential Compton scattering cross section shows a strong density dependence. 
Within a Local Density Approximation, the density dependence of the A propagator results 
in a much better description of the data, as compared to a calculation using the vacuum A 
propagator. Mean-field models which incorporate a reduced effective mass of the nucleon 
and A tend to improve on the results in a Fermi-gas calculation. Both mean-field results are 
quite close indicating that the cross section is not that sensitive to the used effective masses. 
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The present one-body mechanism is unable to describe the data at backward angles. In 
order to improve this, it is imperative that multiple scattering should be incorporated in the 
model. 
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APPENDIX A: ANALYTICAL EXPRESSIONS FOR THE PION SELF-ENERGY 



ReU° ph (q) = 8 



NN 



The real and imaginary part of the pion self-energy can be expressed as 

2 

) [\q\ 2 C (\q\) - C 2 (\q\) + 2M 2 N q 2 Re[L (q)} 
ImU° ph ( q ) = 8 ( ) [2M 2 N q 2 Im[L (q)] 



7, 



where the expressions for the functions Co and C 2 read as 
dk 6(\k\ - k F )9{k F - \k + q\) 



Co(kl) = / 



(2tt) 3 
1 



4E^E~ 

^k^k+q 



E k E k+q 



9Qn 2 \q\ 



(Al) 
(A2) 



(A3) 



3(kT - K - \q] \k\) - 2M 2 + \k\ 2 \ E % + [\q\ 2 - 2(M 2 N + \k\ 2 ) 

k\=k2 

+ \q)\k\ E,*_ ]3] +3E F \k\ 2 + 3\q\M 2 N ln ' 



(E- k + \k\) 



(%H«I + I*|-|$I). 



(A4) 



|fc|=fci 



a 



dk 9(\k\ - k F )9(k F - \k + q\) 



(2tt) 3 
1 

167T 2 |g1 



^k^k+q 



E k ~ E k+q 



E k + E k+q 



(2M 2 N - k 2 F )^\k\ 2 + ^\k\' + ^ [l%T - 40A4|g1 2 - 8M* 



154 - 20M^k 2 F - 30\q\ s \k\ + 10k 2 F \k\ 2 + 20|gl 2 |A;| i + AM 2 N \k\ 2 - 3 



(A5) 



E k 



1 

30 



A\k\ A - \qf + 4il4 + A\q\\k\ 3 - 6\q\ 2 \k\ 2 - \q\ 3 \k\ + SM 2 N \k\ 2 

(%| + 



+ 18M 2 N \q\ 2 + AM 2 N \q\\k\ 



E \k\-\q\ + 



E, 



k\~\q\ 



\k\-\q\ 



\k\=k 2 



|fe|=fci 



(A6) 
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with the integration boundaries 

k\ = max (fop, \q\ — k F ) , k 2 = k F + |g|. 
The relativistic equivalent of the Lindhard function is 



(A7) 



L (q) 



dk 6(\k\ - k F )6(k F - \k + q\) 



L E k- E k+q-<l°- ie 



l€ 

(A8) 



Based on the work in the real and imaginary part of L can be written asQ 



Re[L (q)} 



K2vr) 2 |g1 

- \0{q 2 - AM 2 N ) + d{~q 2 



2E F l{-\q ) + q I { S\q ) + 2 (E^ - E { S ] 

W { - ] U) 



+ 2V^a 



6{AM 2 N - q 2 



arctg 



E 



(+) 



-a 



Im[L (q)] 

where 

I ( +\x) = In 



71 



arctg 



[{2E F + g ) - /(go)] %o - max(0, E^)) 9(E ( + -> - g 



E 



(+) 



-a 



8(27r)2|gl 

- 6{2k F - \q\) [{2E F - g ) - f(q )] 6(q ) 0(-E^ - g ) 



{E ( f> + x)(£f ; - x 



(±) 



(E^ ] +x)(Ei - r 



= E kp+ \j\ ± Ej, 



(±) 



X 



{Ef> + x){E™ - x) 



(±) 



/0»0 = l$l< 



1 + 



4M2 



fc F -|gl 
2 



a = / (?o) 



A, 



|fc|±ki 



£f±2|A;||g1 + |g| 2 + Mjy-. 



(A9) 

(A10) 

(All) 
(A12) 
(A13) 
(A14) 



All analytic expressions have also been checked numerically. 

APPENDIX B: KINEMATICS 

We consider kinematics in the laboratory frame for the 7 A scattering, where the initial 
nucleus is at rest, 



u + M A = J M 2 A + <f 



(Bl) 



1 In the original article |39| some typing errors occur in these formulas. 
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and (u, q) = (k — k',k — k') is the 4-momentum transferred to the nucleus A. The energy of 
the final photon is 

h 

l + fc/M A (l-cos^ 7 ) 1 ' 

where 6 1 is the photon scattering angle. 

The nucleon "effective" momentum p can be found by assuming energy-momentum con- 
servation on the constituent nucleon 



u + y Mjy + p> = \J M% + (p + <f) 2 (B3) 

with the same 4-momentum transfer. 

In the non-relativistic approximation, where |g| <C m, one obtains from Eqs. and 

2pq = f{^f-l). (B4) 

Ma 

A possible solution of Eq. (|B4] ) is given by 

p = —aq (B5) 

with 

„ = !(!_«»)„ 1(1 -I). (B6) 



The component of the momentum p perpendicular to q is not determined from Eq. ( |B4[ ) and 
may be conveniently chosen equal to zero. 

In the relativistic case the solution of Eq. ([B3|) is more complicated. If we seek for an 
effective momentum in the form of Eq. (|B5|) then we obtain two solutions 



^-\^vJ i+2 ^y <B7) 

It is easy to check that for the non-relativistic kinematics ct_ in Eq. ( [B7| ) reduces to Eq. ( |B6| ) . 
Forward scattering (9 1 = 0) is a special case for which u = q = 0. In this case the exact 
solution is given by Eq. flBEp. 

To make a transformation of the single-nucleon T-matrix from the system where the 
nucleon moves with a momentum p to the 7 iV center-of-mass system we consider the 
invariant Mandelstam variables 

s N = (k + pf = M 2 N + 2kE N (p) - 2kp, 

t N = (k - k'f = -2kk\\ - cos# 7 ) . (B8) 

The corresponding center-of-mass 3-momentum and photon scattering angle can be obtained 
from 

\f\ = \k c \ = (s N -M 2 N )/2^, cos9 c N = l + t N /2(k c ) 2 . (B9) 
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FIGURES 




FIG. 1. (a) In the left panel the Pauli-corrected A decay width is depicted as a function of the 
density p (in units of the equilibrium density po [k F = 1.333/m -1 ]) and the A 3-momentum \p&\ 
at IFa = 1232 MeV calculated in the Fermi-gas model, (b) The right panel shows the results for 
this calculation in function of the A 3-momentum \p&\ for the 4 different densities stated in the 
figure. 
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FIG. 2. The imaginary part of the pion self-energy in function of the pion 3- momentum \q\ 
at q° = [(pa) - E F ]/2 for the A kinematics Wa = 1232 MeV and |p A | = 200 MeV within the 
Fermi-gas model and at q° = [(p%)° — Ep]/2 for the A kinematics = M A and |pa| = 200 MeV 
for the mean-field models I and II. The markers on the figures indicate the |q|-values for cosO equal 
to -1 (first marker), (second) and 1 (third) for the specific kinematics stated above which appear 
in the spreading width calculation, see Eq. fl23|). 
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FIG. 3. (a) In the left panel the spreading width is depicted as a function of the density p (in 
units of the equilibrium density po [k° F = 1.333/m -1 ]) and the A 3-momentum |p^| at Wa = 1232 
MeV calculated in the Fermi-gas model, (b) The right panel shows the results for this calculation 
in function of the A 3-momentum \p&\ for the 4 different densities stated in the figure. 
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FIG. 4. The effective nucleon mass Mjy as a function of the density p (in units of the equilibrium 
density po (value depending on the model)) for the mean- field calculations I and II as explained in 



section II D 
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mean-field calculation I 



mean-field calculation II 




FIG. 5. The Pauli-corrected decay width and the spreading width as a function of the 
A 3- momentum \p&\ calculated in the ato— model at = M^(p) for 4 different densities and for 
the two mean- field calculations as explained in the text. 
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FIG. 6. The Pauli-corrected A decay width and the spreading width in function of the 
density p (in units of the equilibrium density po (value depending on the model)) and the A 3- mo- 
mentum |pa| calculated in the mean- field model I and II (marked by index I and II respectively) 
at W£ = M* A (p). 
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7. The calculated cross section for Compton scattering off the proton as a function of 
fixed photon energy, and as a function of photon energy at fixed angle. Data are taken 
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Mean-field 1 




FIG. 8. Differential cross section and photon asymmetry for Compton scattering off 4 He at an 
energy of 206 MeV as a function of angle and at an angle of 37° as a function of energy. The 
curves represent the results of the Local Density Approximation calculation (LDA) and those for 
different densities (in units of saturation density) for the mean-field calculation I. Data are taken 
from |||. 
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FIG. 9. Values of the parameters which define the self-energy of the A resonance evaluated at 
the A invariant mass (Wa) and three momentum (|pa|) appropriate for Compton scattering off 
4 He as shown in Fig. (||). 
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FIG. 10. LDA calculation of coherent Compton scattering off 4 He in the Fermi-gas model and 
the mean-field models I and II. The data are from ref. P,|36|. 
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FIG. 11. LDA calculation of coherent Compton scattering off C in the Fermi-gas model and 
the mean-field models I and II. The data are from ref. |37[| . 
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